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ABSTRACT 

We use the spatial information of our previously published VLT/F0RS2 absorption line spectroscopy 
to measure mean stellar velocity and velocity dispersion profiles of 25 field early-type galaxies at a 
median redshift z ~ 0.97 (full range 0.6 < z < 1.2). This provides the first detailed study of early-type 
galaxy rotation at these redshifts. From surface brightness profiles from HST imaging we calculate 
two-integral oblate axisymmetric Jeans equation models for the observed kinematics. Fits to the data 
yield for each galaxy the degree of rotational support and the mass-to-light ratio M/Ljcans- SO and 
Sa galaxies are generally rotationally supported, whereas elliptical galaxies rotate less rapidly or not 
at all. Down to Mb = —19.5 (corrected for luminosity evolution), we find no evidence for evolution 
in the fraction of rotating early-type (E-I-SO) galaxies between z ^ 1 (63% ± 11%) and the present 
(61% ± 5%). We interpret this as evidence for little or no change in the field SO fraction with redshift. 
We compare M/Ljoans with Af/ivir inferred from the virial theorem and globally averaged quantities, 
and assuming homologous evolution. There is good agreement for nonrotating (mostly E) galaxies. 
However, for rotationally supported galaxies (mostly SO) M/Ljoans is on average ~ 40% higher than 
Af/Lvii-. We discuss possible explanations and the implications for the evolution of M/L between 
z — I and the present and its dependence on mass. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: kinematics and dynamics — galaxies: 
evolution — galaxies: fundamental parameters 



1. INTRODUCTION 

Detailed studies of nearby early-type galaxies pro- 
vide insight into their stellar populations, dark mat- 
ter content and kinematic and spatial structure. In 
particular, resolved kinematic dat a reveal the under- 
lying gravitational po tential (e.g.. Ivan der Mareilll991l : 
iCappellari et al.l I2006D. the presence of super- massive 
black holes (e.g.. iRichstone et al.1 llQQOl: [V erolm e et al.l 
l2002f l and dark halos (e.g., Kronawit ter et al.ll2000f l. and 
the relative contributi ons of p ressu re and rotation to 
its orbital energy (e.g ., lBinnevlfl978l : iDavies et al.lll983l : 
lEmsellem et aljfcoOTh . 

Recently, the first, probing steps have been made 
toward spatially resolving the stellar motions in more 
distant early-type galaxies using absorption-line spec- 
troscopy. The past few years have seen stud- 
ies aimed a t constraining the da rk halos of lens- 
ing galaxies (jKoopmans et al.l l2006f ). dyn amically dis- 
tingu ishing cluster E and SO galaxies (|Moran et al.l 
l2007tl. con strain ing the evolution in rotation rate 
(|van der Marel fc van DokkumI l2007af) . and constrain- 
ing the evolution of the mass-to-light ratio {M/L) us- 
ing the M/L versus a relation at inte rmediate redshift 
(|van der Marel fc van DokkumI l2007bf ). Obviously, the 
spatial resolution and errors in the kinematic quantities 
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Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS 5- 
26555. 
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are much larger than for local galaxies, and consequently, 
the level of detail is far lower. Still, the spatially resolved 
information that is obtained by going beyond studies of 
global quantities alone, combined with the addition of 
cosmic time to the equation, has already revealed sev- 
eral interesting results. 

l^Ioran et al. (2003) used rotation curves to address 
the question whether E and SO galaxies are coeval or 
not. They found that the relative number of rotation- 
ally supported early-type galaxies at z ~ 0.5 is lower 
than today, especially in clusters. Combined with ex- 
tensive studie s of vi sually classifie d samples of galaxies 
out to z ^ 1 (iDressler et al.lll99"7t iPostman et al.l 120051 : 
ISmith et al.ll2005l ). this is the best evidence to date for 
relatively recent SO formation, presumably through the 
transformation of in-falling, star-forming Sa-like galax- 
ies into quiescent cluster SO galaxies (e.g.. lGunn fc Got! 
[T97llLarson et al.lll98l . 

This classical picture does not address the high fraction 
of SO galaxies at low densities. Approximately 60% of all 
early-type galaxies in the field are SO galaxies. Moreover, 
recently it was shown that the total fraction of early- type 
galaxies does not significantly ch ange between z ^ 0.8 
and the pre sent either in clusters ([Holden et al.ll2007l l or 
in the field ()van der Wei et al]|2007f ) if only galaxies with 
masses M > 0.5M* are considered. This calls into ques- 
tion the claim that SO galaxies evolve differently from E 
galaxies, although we note that this has so far not been 
explicitly addressed. Because of these reasons it is im- 
portant to not lose sight of alternative scenarios for the 
formation of SO galaxies. Merging that does not com- 
pletel y destroy th e disk is another proposed mechanism 
(e.g., Bckki 1997). This is more similar to the favored 
formation mechanism for E galaxies and would therefore 
imply that E and SO are in fact one class of objects with 
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a large range in bulge-to-disk ratios. This ratio may be 
related t o the mass ratio and the dissipation of the merg- 
ers (e.g.. iNaab et al.ll2006al ). In this scenario E and SO 
galaxies are expected to be roughly coeval. 

A nother interesting observation i n thes e contexts is 
that Ivan der Marel fc van DokkumI ([2007a) found some 
evidence that the relative number of rotationally sup- 
ported cluster elliptical galaxies decreases slightly be- 
tween z ~ 0.5 and the present. This points to- 
ward a scenario in which elliptical galaxies gradu- 
ally lose angul ar momentum throu gh interactions and 
mergers (e.g., iNaab et aD l2006bl ). This result is 
not necessar il y inc onsistent with the conclusions from 
iMoran et al.l ()2007[ ). The latter authors used rotation 
to trace the number of SO galaxi e s, w hereas the sample 
of Ivan der Marel fc van DokkumI (|2007a) contained (al- 
most) no SO galaxies by construction. Therefore, these 
studies probed different aspects of the early-type (E-I-SO) 
galaxy population. 

Besides constraining rotation, resolved kinematic data 
also allow for a measurement of M/L. The evolution of 
M/L with redshift is typically studied using estimates 
based on the virial theorem or the funda mental plane 
(e.g.. Ivan Dokkum fc van der Marel I2007D . This uses 
only globally averaged quantities (luminosity, radius, in- 
tegrated velocity dispersion), combined with assump- 
tions about the dynamical structure of the galaxies and 
the homology of their evolution. Spatially resolved kine- 
matics have the advantage that they allow the M/L to be 
inferred through detailed dynamical modeling, which re- 
moves some of the assumptions inherent in these studies. 
Systematic effects may affect not only the measured, av- 
erage evolution of M/L but also the evolution of the tilt 
of the fundamental plane, i.e., the dependence oi M/L 
evolution o n galaxy mass (e.g.. Ivan der Wei et al.|[2005l : 
iTreu et al. 2005). Van der Marel fc van Dokkum (2007b) 
addressed this for a sample of clusters elliptical galaxies 
at z '^ 0.5. They found that bright, nonrotating galaxies 
are consistent with homologous evolution, but that the 
situation may be more complicated for less luminous, ro- 
tating galaxies. However, the small size of their sample 
and the limited number of SO galaxies left room for dif- 
ferent interpretations of this result. 

In this paper we study the resolved kinematics of 
a sample of field early- type galaxies at z ~ 1, using 
spatially resolved spectra from the Very Large Tele- 
scope (VLT) and imaging from the Hubble Space Tele- 
scope (HST). Th i s dat aset was presented earlier by 
Ivan der Wei et all ()2005l ). Our analysis methodology is 
largely similar to that of van der Marel fc van Dokkum 
(2007a, 2007b). However, our sample is different from 
theirs in three important areas. First, we focus on field 
galaxies instead of cluster galaxies. Second, our sample 
includes significant numbers of SO and Sa galaxies, in- 
stead of just E galaxies. And third, we work at z ~ 1 
instead of z ~ 0.5, thus providing a much extended base- 
line for studying evolution. In particular, we quantify 
the evolution of the relative number of rotationally sup- 
ported field early- type galaxies from z ^ 1 to the present, 
and we examine the difference between the M/L inferred 
from spatially resolved and unresolved kinematic data. 

We have organized the paper as follows. In Sj2|we in- 
troduce the data. In fj3| we briefly describe the dynam- 
ical modeling. In Sj4| we present our results regarding 



rotational support. In El5|we present our results regard- 
ing the M/L. In f}6| we summarize the main conclu- 
sions. We adopt the following cosmological parameters: 
(l^Af, ^A, h) = (0.3, 0.7, 0.7). 

2. DATA 

2.1. Spectroscopy 

Very deep VLT/F0RS2 spectroscopy of a sample of 
magnitude-limited, morphologically selected field early- 
type galaxies at redshifts 0.6 < z < 1.2 in the Chan- 
dra Deep Field-South (CDF-S) and the foreground of the 
z = 1.24 cluster RDCS 1252 9-292 7 (CL1252) were pre- 
sented bv lvan der Wei et al.l (|2005D . That paper provides 
a full description of the data reduction, including the 
technique used to measure the global velocity dispersion 
of each galaxy from its absorption features. In this pa- 
per we examine a subsample of 25 early-type galaxies (10 
E, 9 SO, and 6 Sa galaxies) that have such high-quality 
spectra that variations of mean velocity v and velocity 
dispersion a along the slit can be measured. The median 
redshift of this subsample is z = 0.97. Since this is a 
magnitude-limited subsample it is not biased against ro- 
tating or nonrotating galaxies, and it is representative for 
the entire population of E-I-SO galaxies at this redshift. 

For high-redshift galaxies, a is usually determined from 
an extracted spectrum that is averaged over multiple pix- 
els in the spatial direction in order to optimize for signal- 
to-noise ratio (S/N). Here, instead, we measure v and 
a for each individual row of pixels (0.251ong) for which 
the S/N is sufficiently high for a robust measurement. 
We use t he same wavelengt h range and template stars as 
used by Ivan der Wei et all (|2005D such that differences 
between the results from averaged and resolved spectra 
are solely due to differences in the observed spectra them- 
selves. The mean velocity profiles, typically measured for 
5-7 pixel rows, are shown in Figure [T] It is readily appar- 
ent that many galaxies show signs of significant rotation. 
Measuring a requires better S/N, therefore a can usu- 
ally only be measured for 3-5 pixel rows. The S/N is 

typically 15 A for the outermost pixel rows for which 
we measure a (the central pixel rows obviously have the 
highest S/N). For two galaxies (CDFS-18 and CDFS-23) 

we have S/N < 8 A only for the outermost pixel rows, 
which is not ideal for velocity dispersion measurements. 
However, for consistency we retained these rows in our 
analysis because the large resulting error bars are explic- 
itly accounted for in our modeling described in fj3l In 
Figure m we show the a profiles. In both Figures [T] and |2] 
we over-plot the best-fitting models that we describe in 
Sj3| below. 

2.2. Photometry 

GOODS^ provides deep, publ icly available HS T /ACS 
imaging of the CDF-S dCiavalisco et al.ll200l . ACS 
imaging for the fiel d of CL1252 is also available 
(jBlakeslee et al.l [20031 ). We used these data to assign 
morphological classifications (E, SO, Sa) based on visual 
inspe ction of the i mages, following the strategy outlined 
bv P ostman et all ([2005). These classifications are used 
throughout this paper. 

^ http://www.stsci.edu/scicnce/goods/ 
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Fig. 1. — Profiles of mean velocity v for 25 field early- type galax- 
ies at redshifts 0.6 < z < 1.2. The points indicate the measured 
velocity at each spatial pixel in the spectroscopic slits. The solid 
curves show the predictions of the best-fitting models. For compar- 
ison, dotted curves are the corresponding predictions when seeing 
convolution and pixel/slit binning are not taken into account. The 
visual morphology, the redshift, and the inferred rotation param- 
eter k are given for each object. An asterisk indicates that the 
maximum value for k that is physically possible is adopted instead 
of the best-fitting value for k. The ID numbers correspond to those 
given in van der Wei ct al. (2005), where more information regard- 
ing the photometric and kinematic properties can be found. 

Van der Wei et al. (2005) used the F850LP images to 
fit de Vaucouleurs profiles to the objects in their sam- 
ple in order to determine their effective radii R^s and 
surface brightnesses /Lteff- For the present analysis we 
use the ACS data to measure the full surface brightness 
profile for eac h galaxy. We follow the same procedure de- 
scribed by va n der Marel fc van Dokk um (2007a). First, 
the images are deconvolved with the point-spread func- 
tion (PSF; for which we use stars in the field) using 
the Lucy-Robertson algorithm. Subsequently, we fitted 
elliptical isophotes to the two-dimensional images with 
the IRAF task ELLIPSE. This yields one-dimensional 
profiles of major-axis surface brightness, position angle, 
and ellipticity. For galaxies at redshifts z > 0.85 we 
use the ACS F850LP images, and for galaxies at red- 
shifts z < 0.85 we use the ACS F775W images. With 
these choices, the observed wavelength corresponds as 
closely as possible to the rest-frame B band. The ob- 
served surface brightness profiles are transformed into 
the rest- frame B band using t he m ethod presented in 
Ivan Dokkum fc van der Marell (|2007f) . This facilitates 
the comparison with local galaxy samples. 

3. MODELS 

For a full description of the modeling proce dure we re- 
fer to Ivan der Marel fc van DokkumI (|2007a[ ). In short, 
the rest-frame i?-band photometry described in ii2.2l is 
fitted with the projection of a parameterized, oblate ax- 
isymmetric, constant axial-ratio luminosity distribution. 
The results of the modeling are fairly insensitive to vari- 
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Fig. 2. — Profiles of velocity dispersion a for the same 25 field 
early-type galaxies at redshifts 0.6 < z < 1.2 as in Fig. [T] The 
points indicate the measured dispersion at each spatial pixel in 
the spectroscopic slits. The solid curves show the predictions of 
the best-fitting models. For comparison, dotted curves are the 
corresponding predictions when seeing convolution and pixel/slit 
binning are not taken into account. The visual morphology, the 
redshift, and the rest-frame S-band model A//L in solar units are 
given for each object. 

ations in the unknown inclination (van der Marel fc van 
Dokkum 2007a, 2007b). In the discussion below we adopt 
for each galaxy the inclination angle that is most likely, 
given the probability distribution of intrinsic axial ratios 
derived from large galaxy catalogs of the local universe. 
Given the three-dimensional luminosity density, the 
Jeans equations are solved under the assumptions of a 
constant M/ L and a two-integral distribution function 
/ = f(E, L^), where E is the energy and L^ the angular 
momentum around the symmetry axis. The models have 
v'j^ = v1, so their velocity distribution is isotropic in the 
meridional plane. The ratio of i;? to v\ is determined 
by the requirement of hydrostatic equilibrium. The sec- 
ond azimuthal velocity moment is split into mean and 
random components according to the convenient param- 
eterization 
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A value fc = yields a galaxy that is nonrotating and 
fully pressure supported. A value |fc| = 1 yields an oblate 
isotropic rotator, i.e., a galaxy with sufficient rotation to 
fully explain its flattening. The value of |fc| can be larger 
than unity, with a physical maximum that depends on 
the shape of the galaxy, and set by the requirement that 
cr? = u? — {v^Y is everywhere positive. In the following 
we always choose the major axis position angle (which is 
indeterminate module 180°) so that k is positive. 

The solutions of the Jeans equations are integrated 
along the line of sight, and then convolved in luminosity- 
weighted sense with the seeing (typically 0.6-0.8) and 
with the size of the pixels (0.251ong) in the slit (l.Owide) 
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Fig. 3. — B-band luminosity vs. rotation rate. The large sym- 
bols are the 25 field early-type galaxies at redshifts 0.6 < z < 1.2, 
for which we show the rotation parameter k. The luminosities are 
corrected for 0.58 dex of evolution b etween z = 1 and the present 
l|van Dokkum &: van der Marelll2007I V Filled symbols are E galax- 
ies, open circles are SO galaxies, and open, dotted circles without 
error bars are Sa galaxies, as determined by visual classification. 
The local comparison sample is shown as small symbols, with the 
different morphologies distinguished by the same symbols as used 
for the distant sample. (V/cr)*, which is comparable to fc, is used 
to quantify the rotation rate for the local galaxy sample. The dis- 
tributions of the local and distant samples are similar. 

through which the galaxies' spectra were taken, taking 
the angle between the orientation of the slit and the 
major axis into account. The match of the predictions 
to the kinematic data from ^ 32. II is optimized by finding 
the best-fitting value of the two available model param- 
eters. The mass-to-light ratio AI/L determines the total 
amount of rms motion (vV + ct^) whereas the parame- 
ter k determines the amount of rotational support (which 
is related to v/a). Solid curves in Figures [1] and [3] show 
the best model fits. For comparison, dotted curves show 
the model predictions if no account is taken of seeing and 
pixel/slit binning. Clearly, the latter makes a consider- 
able difference and it is therefore important that this is 
properly modeled. Best-fit values of k are listed in the 
individual panels of Figure [1] 

4. ROTATION AT Z ~ 1 
4.1. Evolution in Rotation Rate 

As is readily apparent from Figure [1] many early- type 
galaxies at 2 ~ 1 show signs of rotation. In Figure [3] 
we show using large symbols the rotation parameter k as 
a function of absolute _B-band magnitude, distinguishing 
between galaxies with different visual morphologies. The 
absolute magnitudes are corrected for luminosity evolu- 
tion u sing the results of Ivan Dokkum &: van der Marell 
()2007l ). Slightly different values would have been rea- 
sonable as well (see, e.g., ^5^, but none of our results 
depend sensitively on this. 

Only five out of the 25 galaxies show little or no ro- 
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Fig. 4. — EUipticity vs. rotation rate [k for the distant sample, 
(V/cr)* for the local sample]. The symbols that indicate different 
visual classifications are the same as in Fig. |3] The ellipticity dis- 
tribution, like the rotation rate distribution, is similar for the local 
and distant samples. 

tation. These are all elliptical galaxies. By contrast, all 
SO and Sa galaxies show significant rotation, with many 
consistent with being rotationally supported. Figure [5] 
shows the rotation parameter /c as a function of the ap- 
parent ellipticity e = 1 — (b/a) , using similar symbols as 
in Figure [3] This shows that the nonrotating galaxies 
tend to be rou nder than mo s t of t he rotating galaxies. 

As noted bv lMoran et al.l (|2007f) . a single-component, 
spheroidal model may not be representative for SO 
galaxies, some of which clearly show disks. Our 
models can in principle b e extended to include disks 
(jCinzano fc van der Marelll9 94). but we have not ex- 
plored that here. Therefore, the inferred k may not al- 
ways necessarily correspond to the physical property it 
represents in the model. Nonetheless, as a fit parameter, 
k still provides an effective measure of the relative impor- 
tance of rotation in the galaxy. In fact. Figures [3] and [5] 
show that E and SO galaxies can be distinguished fairly 
successfully based on their kinematics, as quantified by 
the parameter k. Additional use of luminosity and ax- 
ial ratio information can further increase the accuracy of 
such a kincmatical classification scheme, as advocated by 
iMoran et al.l (|2007h . The only galaxy that seems some- 
what out of place in Figure [3] is the most rapidly rotating 
elliptical galaxy. However, this is also the least luminous 
galaxy in the sample, and it may well have been visually 
misclassified. 

The clean kinematic separation between different mor- 
phological types at 2 ~ 1 is surprising gi ven the complex 
situation for local early-type galaxies (jEmsellem et al.l 
l2007f ). Both E and SO galaxies are of ten fast rotators 
with overlapping kinematic properties (jCappellari et all 
|2007| ). The clean separation at higher redshifts is there- 
fore quite likely partially artificial. Morphological clas- 
sifications are difficult at high redshift, and recent work 
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Fig. 5. — Histograms of early-type (E+SO) galaxy rotation rates. 
The parameter k is used for the distant sample (solid line with 
Poisson error bars) and the quantity (V/cr)* is used for the local 
sample (dashed line). Only galaxies brighter than Mg = —19.5 
are included (corrected for luminosity evolution, in the case of the 
distant sample). There is no significant difference between the two 
samples in terms of relative numbers of slow-rotating and fast- 
rotating galaxies. 

has shown that the relative number of SO galaxies at high 
redshift is systematically underestimated as evidenced by 
the lack of SO galaxies with e < 0.3 (B. P. Holden et al, in 
preparation). In addition, elongated, rotating elliptical 
galaxies can be misclassified as SO galaxies. Nonetheless, 
these concerns are not relevant for the main goal of this 
paper, the evolution of the rotation of early-type galaxies 
as a single class of objects. 

To address the issue of evolution in the rotation rate, 
we compare the z ~ 1 sample with a local sample ex- 
tracted from the L yon Extragalactic DataBase (LEDA; 
iPaturel et al.lll997l ) , which is distribute d and made avail- 
able as HyperLeda (jPaturel et al.ll2003[ )^. We restrict the 
local sample to galaxies within ^ 40 Mpc (distance mod- 
ulus m — M < 33). In order to construct a field/group 
E+SO sample we included only galaxies with fewer than 
80 neighbors brighter than Mb — —19.5 within a cylin- 
der with a 5 Mpc diameter (at the distance of the galaxy) 
and with a difference in radial velocity with respect to the 
cosmic microwave background of less than 1000 km s~^. 
This effectively selects a complete sample of 179 galax- 
ies (70 Es, 109 SOs) outside overdense structures (most 
notably, the Virgo Cluster). The median projected sur- 
face density is E = 0.4 Mpc~^, which is derived from 
the distance to the 7th nearest neighbor down to the 
luminosity limit of the magnitude-limited sample with 
the same radial velocity within 1000 kms~^. This is 
the same value as for a l arge, volume-limited sam ple of 
galaxies at z '-^ 0.03 (see Ivan der Wei et al.l 12007^. It is 
also comparable with the estimated surface density (in 
co-moving units) for a sample of distant field galaxies 

^ http://lcda.univ-lyonl.fr 



(see Ivan der Wei et al] l2007f ) of which our distant sam- 
ple is a sub-sample. Maximum stellar velocities V and 
velocity dispersions g ar e known for 55 Es and 43 SOs 
(jPrugniel fc Simieiilll996D . There are 37 galaxies in this 
sample for which multiple, independent measurements of 
the rotation velocity V are available. For 31 of these the 
measurements are consistent within the errors. The mea- 
surements are inconsistent only for t he remaining six (see 
Table A2 of lPrugniel fc Simienlll996h . Because especially 
galaxies with peculiar properties tend to be targeted mul- 
tiple times, the measurements for the remainder of the 
sample may actually be even better. Therefore, the kine- 
matic measurements for the local sample are generally 
reliable and reproduceable. 

We quantify the rotation rate of the local sample using 
the quantity {V/a)* = (y/a)A/(r^^7)7e. This approxi- 
mates the ratio of the actual rotation rate of the galaxy 
to the ratio expected for an o blate isotro pic rotator of 
the given observed axial ratio (jDavies et a l. 1983). This 
quantity is therefore similar to our model parameter k. 
The only difference is that k is defined locally in the 
meridional plane by equation (1), whereas (V/a)* is de- 
fined in terms of the globally defined projected quantities 
V, a and e. Van der Marel (1991; see his Fig. 3a) made 
a direct comparison of k and (V/a)* for a sample of local 
galaxies and found a good agreement to ~10%. For the 
present paper we therefore compare {V/a)* for the local 
sample directly to k for the distant sample. 

The local sample is shown in Figures |3I and |3I us- 
ing small symbols. As compared to the distant sample, 
there are some differences in the types of galaxies that 
show high rotation rates. In the distant sample these 
are almost all fiattened, visually classified SO galaxies. 
By contrast, in the local sample these include quite a 
few roundish, visually cla ssified E galaxies. It has been 
known for some time fe.g.. lRixfc White 1990) that local 
galaxy catalogs are deficient in such galaxies. Of course, 
distant samples may be plagued by similar, or worse, 
problems. However, our sample size is too small to make 
definite statements about this issue. 

As noted earlier, possible uncertainties in the relative 
classifications of E and SO galaxies do not affect the as- 
sessment of the rotation rate distribution of early-type 
(E-t-SO) galaxies in general. Figure |5| shows that the ro- 
tation rate histograms are statistically indistinguishable 
for the local and distant samples: 63% ± 11% (12 out 
of 19) of the E+SO galaxies in the distant sample have 
high rotation rates (fc > 0.6), and 59% ± 5% (58 out of 
98) of those in the local sample have high rotation rates 
[{V/a)* > 0.6]. 

The similarity between these two numbers is striking, 
however, they may not be directly comparable as the 
kinematic data for the local sample are not complete, 
nor representative for the morphological composition of 
the population: an E galaxy is twice as likely to have 
kinematic data as an SO galaxy (see the values above). 
If we, very crudely, assume that all SO galaxies rotate 
and all E galaxies do not, then the fraction of rotating 
E+SO galaxies would be 69%. However, in reality the rel- 
ative fractions of rotating/nonrotating E and SO galaxies 
are not very different (55% and 65%, respectively, see 
the small data points in Fig. [3]). Based on these relative 
numbers of rotating galaxies for E and SO galaxies as 
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separate classes we estimate that the true relative frac- 
tion of rotating early-type galaxies is 61%, only slightly 
higher than the measured value of 59%. 

Concluding, we find no evidence for evolution in the 
rotation rate of field early-type galaxies between z = 
1 and the present down to the luminosity limit of the 
distant sample, Mb = —19.5. The distant sample is 
small (as evidenced by the large error bars in Fig. [5]) , but 
we would have detected an increase or decrease by more 
than 25% in the relative number of rotating galaxies if 
there were such evolution. 

4.2. Comparison with Previous Results 

The only previous measurement of the evolution of 
the relative number of rotationally supported early-type 
(E-I-SO) galaxies is that of Moran et al. (2007) . Those 
authors found substantial evolution for cluster early-type 
galaxies in the sense that fewer early-type galaxies are ro- 
tationally supported at z ~ 0.5 than locally, which is in- 
terpreted as evidence for a lower SO fraction in the distant 
clusters. This result for cluster galaxies is not directly 
comparable to ours. However, t he field early-type galaxy 
sample from lMoran et all (|2007l ) contains only 43%±14% 
rotating galaxies, which is less than the value for local 
galaxies inferred from the LED A sample described above. 
If this (marginally significant) decrease between z ~ 0.5 
and the present is real, more evolution between z ~ 1 
and the present can be expected, but in our sample we 
do not observe this. Besides the low significance, this 
comparison is furthermore hampered by differences in 
methodology. Therefore it is interes ting to compare our 
results and methods with those from lMoran et all ()2007l ) 
in more detail. 

Both we and iMoran et al.l (|2007l ) measured spatially 
resolved kinematical profiles of gala xies, but beyond that 
our methods are very different. Mor an et al.l ([2007) infer 
global quantities (the rotation v, velocity dispersion a, 
and ellipticity) , and they do so directly from the data. 
Observational effects are not taken into account, and the 
internal dynamics of the galaxies are not modeled. Our 
method aims to model the internal dynamics, to take 
into account known differences between galaxies (e.g., in 
their surface brightness profile) , and to account explicitly 
for all known observational inf luences on th e mea sured 
quantities. The approach of iMoran et al.l ()2007t ) has 
the advantage of being simple and model- independent. 
However, ignoring the effects of seeing convolution and 
pixel/slit binning can bias the results. Figure [U shows 
that these effects have a significant impact on both the 
measured rotation value at the outermost radius and the 
maximum value of the rotation curve. As a result, the 
relative number of rotating galaxies can be underesti- 
mated. 

To quantitatively illustrate the effect of seeing convo- 
lution and pixel/s lit binning, w e apply one of the kine- 
matic classifiers of IMoran et al.l |2007), namely v/{l — e), 
to our datasets. Here, v is defined as half the velocity 
range of a fitted straight line from end to end of the 
measured velocity profile, approximating the maximum 
rotation velocity. From the LEDA sample we derive that 
65% ± 5% of the local p opulation satisfies the criterion 
from lMoran et all (|2007l ). virtually the same fraction as 
galaxies with (V/a)* > 0.6 (see H4.ip . This demonstrates 
that for the purpose of distinguishing rotating and non- 



rotating galaxies, the two methods are in principle com- 
parable. 

In our z ~ 1 sample only 7 out of 19 E-I-SO galax- 
ies (37% ± 11%) ha ve v/(l - e) > 85 km s~\ consistent 
with the results from M oran et all (|2007l ) for z ~ 0.5 field 
galaxies, but substantially less than the number of galax- 
ies with k > 0.6 (63% ± 11%) in our z ~ 1 field sample. 
However, if we use the maximum model line-of-sight ve- 
locity, with observational effects taken into account (i.e., 
the dotted curves instead of the solid curves in Fig. [I] 
then we find a higher fraction with v/{l — e) > 85kms^^: 
63% ± 11%, the same as the fraction of galaxies with 
k > 0.6. This implies that omission of the effects of see- 
ing convolution and pixel/slit binning can lead to signif- 
icant underestimates of the number of rotating galaxies 
at high redshift. This may well explain the relatively low 
nu mber of kinematical ly classified field SO galaxies found 
bv IMoran et al.l ()2007D at z '^ 0.5. However, we cannot 
assess the question as to what extent the above-described 
bias contributes to the evolution that lMoran et al.) (J2007T ) 
found for cluster galaxies. 

Interesting ly, th e kinematic classifications by 
IMoran et al.l (|2007D agree, statistically speaking, 
with visual morphological classifications, at both low 
and high redshifts. This argues either against a bias 
in their results or that visual classifications can suffer 
from a similar bias against rotating/SO galaxies. An 
indication that high-z classifications of SO galaxies are 
indeed hampered by systematic problems is that the 
ellipticity distribution for the cluster early-type galaxy 
population does not change significantly with redshift, 
suggesting that high-z SO galaxies are systematically 
misclassified as E galaxies, more so than the other way 
around (B. P. Holden et al., in preparation). Still, to 
what extent biases contribute to an apparent decline 
in the SO fraction with redshift will probably continue 
to be a matter of debate for some time to come. Also, 
the answer may well be different for cluster and field 
galaxies. However, it is certainly intriguing that our 
detailed observations and modeling of field galaxies at 
z ~ 1 show no evolution in the fraction of rotating 
early- type galaxies. 

Note that our result does not necessarily imply that 
there is no evolution in the early-type galaxy popula- 
tion. On the one hand, rotationally supported early- 
type galaxies may merge and produce pressu re-supported 
elliptical galaxies (e.g.. JNaab et al.ll2006bl ). This pro- 
cess would produce an increase in the fraction of elon- 
gated/rotating galaxies with redshift, for which some 
tentative evidence exists in cl u sters fr om the results of 
Ivan der Marel fc van DokkumI (j2007aD . On the other 
hand, rotating, quiescent galaxies may be formed out 
of Sa galaxies that cease to form stars, contributing to 
the i ncrease in the st e llar mass density o f red galaxies 
(e.g., iBell et al.l 12004 iBrown et all |2007| ) and the de- 
crease of the cosmic average star form ation rate (e.g., 
iLe Floc'h et al.! '2005: Nocske et al.'20Q3). As long as the 
mechanisms that increase and decrease the relative num- 
bers of rotating and nonrotating galaxies do not change 
their ratio, continued growth of both the E and the SO 
population is possible up to the present day. 

Neither the resu lts of our own study nor that of 
iMoran et al.l (I2007D are directly comparab l e to th ose ob- 
tained bv Ivan der Marel fc van DokkumI ()2007al ) using 
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the same technique as was in this paper. This is be- 
cause their study at z ^ 0.5 dealt with a cluster sample 
that was preselected to contain almost no SO galaxies. 
Therefore, it can shed little light on the evolution of the 
SO fraction with redshift. 

5. MASS-TO-LIGHT RATIOS 

5.1. Comparison between Model and Virial M/L 
Estimates 

According to the virial theorem, the mass of a stellar 
system can be written as 



1 - 



Mvir = 



G ' 



(2) 



where RoS is the effective radius and Coff is 
the luminosity-weighted velocity dispersion measured 
through an aperture of size i?cff • The homology parame- 
ter (3 generally depends on the detailed density structure 
and velocity dispersion anisotropy of the galaxy. Plausi- 
ble models can span a wide range of /? values. However, 
ICappellari et al.l (|2006f ) found that early-type galaxies in 
the local universe all follow /3 = 5.0 ± 0.1, with surpris- 
ingly low scatter. This calibration was obtained by calcu- 
lating galaxy masses M from detailed dynamical models 
for integral-field kinematical data, including higher-order 
velocity moments. The quantity aes was measured by 
integration over the spatially resolved two-dimensional 
velocity field. The virial mass-to-light ratio M/Lvir can 
be calculated upon division by the luminosity L. The 
latter can be measured through aperture photometry or 
it can be estimated as L = 27ri?^ff/off, where IcS is the 
average intensity inside -Roff • 

Equation ^ is commonly used to estimate the mass of 
distant galaxies. In doing so, one generally measures the 
velocity dispersion cTc through a spectroscopic aperture. 
The velocity dispersion Ceff is then estimated by apply- 
ing a correction formul a that is based on typi cal observa- 
tions of local galaxies (iJorgensen et al.lll995D . However, 
the findings of 'Ca ppellari et all ( 2006D do not guarantee 
that equation ^ with /3 = 5 is as accurate as it is in the 
local universe, for at least two different reasons. First, 
galaxies may evolve so that the homology parameter /3 
could be a function of redshift. Second, the accuracy 
and applicability of the corrections from Cc to Ucff are 
not guaranteed at high redshift. The spectroscopic aper- 
tures and the seeing in use there often exceed the galaxy 
size i?cff- Observations that adequately mimic this are 
not available in the local universe. Moreover, even in 
the local universe croff/cc shows strong variations from 
galaxy to galaxy fCappellari et al. 2006). It is therefore 
necessary to calibrate equation ^, and the appropriate 
value of /3, at high redshift in the same way as was done 
locally. This is possible for our sample by comparing the 
M/Ljeans from our dynamical models {^j with the val- 
ues of M/Lvir- Our data quality is obviously not compa- 
rable to the two-dimensional velocity fields and higher- 
order moments that are available locally. Nonetheless, 
this calibration is unique and has not previously been 
done at these redshifts. 

In Figure [S] we compare Af/Lvir and M/L jeans for our 
sample at z ~ 1. We compute d M/Lyjr using P = 5, with 
(7c and Rcs as measured by Ivan der Wei et al.l |2005). 
The spectroscopic aperture was defined by the width 
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Fig. 6. — Jeans M/L vs. virial M/L (calculated using /3 = 5 in 
eq. [2]) for the 2 ~ 1 galaxy sample. The symbols that indicate 
different visual classifications are the same as in Fig. [3] E galaxies 
follow the relation Af/Ljcans ~ Af/Lvir (dotted line). However, 
SO galaxies on average have A//L jeans systematically higher than 
Af/Lvir by ~ 40%. The solid line is the least-squares fit for the 
full sample, which has M/L j^ans oc Af/L^jJ?^ (and has Af/L jeans 
systematically higher than Af/Lvir by ~ 20%). 

of the spectroscopic slit (1") and the height of the ex- 
tracted spectrum (usually, 1.25 "). The ratio aes /o'c was 
estimated using the formula of ICappellari et all (20061. 
which is consistent with those of lJorgensen et all (1995). 
As for local galaxies, the slope of the best-fitting linear 
relation between M/Lvir and M/Ljeans is consistent with 
unity. However, there is a systematic offset that appears 
to depends on the galaxy type. This is seen in the resid- 
uals log[(M/L)joans/(M/L)vir] showu in FigurejTj For E 
galaxies, the residuals are consistent with zero. There- 
fore, /3 = 5 is appropriate for E galaxies at high redshift, 
and as a class E galaxies are consistent with homologous 
evolution. For SO (and Sa) galaxies, on the other hand, 
M/L Joans is on average almost 40% higher than M/Lvir- 
As shown in Figures jTJa and |3), these galaxies tend to 
be flatter and more rapidly rotating than the E galax- 
ies . These results are consistent with the earlier findings 
of Ivan der Marel fc van DokkumI (|2007br ) in clusters at 
z '^ 0.5. However, the trends are much clearer in the 
present sample because of its higher fraction of rapidly 
rotating SO and Sa galaxies. 

There are several possible explanations for the trends 
observed in Figure jT] First, SO and Sa galaxies at z ~ 1 
may have /3 > 5. However, this would imply that these 
galaxies evolve nonhomologously, while E galaxies do 
not. This seems to us unlikely as in that case it would 
be an odd coincidence that /9 = 5 for all types of lo- 
cal galaxies. Second, it is possible that our models have 
overestimated the M/Ljcans for rapidly rotating galaxies. 
Such galaxies may have rapidly rotating disks that are 
not well represent ed by our constant axial ratio o blate 
models. However, ICinzano &: van der Marel (|1994f ) ad- 
dressed this issue explicitly for the well-studied nearby 
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Fig. 7. — The logarithm of the ratio {M/L)jcans/(^/i)vir vs. eUipticity {left panel), v/a with v as defined by IMoran et al.l II2007I ) 
(middle panel), and efli'ective radius (right panel). The symbols that indicate different visual classifications are the same as in Fig. |3] The 
residuals tend to be positive for galaxies of SO or Sa type, for galaxies with significant eUipticity or rotation rate, or for galaxies of small 
size. Possible causes for this, and implications for studies of fundamental plane evolution, are discussed in the text. 



early-type disk galaxy NGC 2974 and found that a single- 
component model overestimated the AI/L by only 12%. 
This is not sufficient to explain the trends seen in Fig- 
ure[71 Third, it is possible that the virial equation ^ for 
M/Lvh- underestimates the true mass of rapidly rotating 
galaxies at high redshift. Rotation contributes to the hy- 
drostatic support of the galaxy, so if it is ignored then 
the M/L will be underestimated. While rotation does in- 
crease CTcff through the effect of line broadening, this may 
not fully capture the true hydrostatic importance of the 
rotation component. This may be especially important if 
the seeing is of order of or larger than the scale at which 
rotation is manifest. Such problems are avoided for local 
galaxies for which much higher quality data are available. 
This may explain the apparent differences between the 
M/L of rotating galaxies at low and high redshifts. For 
our distant galaxies, correction formulae for (Tcs/cTc are 
used that are not themselves calibrated at high redshift. 

5.2. Implications for the Evolution of M/L 

In the previous section we demonstrated that Af/Lvir 
may be too low for rapidly rotating galaxies at higher 
redshifts, but not for local early-type galaxies. Since the 
fundamental plane is in essence a correlation between 
Af/ivir and other global galaxy parameters, this will 
bias the amount of M/L evolution inferred from funda- 
ment al plane studies (see Ivan Dokkum &: van der Marell 
120071 and references therein). Van der Wei (2005) 
found that the rate of luminosity evolution for the sam- 
ple of field galaxies that we also use in this paper is 
Alog(M/L) = (-0.76 ± 0.07)0. Figures [6] and [t] sug- 
gest that, averaged over the rotating and nonrotating 
galaxies in our sample, M/Lvir rnay be too low by 
^ 20%. If so, then the evolution in M/L is reduced 
to log(M/L) = (-0.69 ± 0.07)z. This is steeper than 
the evolution found for cluste r galaxies, but not by much 
(|van Dokkum fc van der Marel.2007 ). 

The difference between M/Ljeans and Af/Lvir depends 



on the galaxy rotation rate. Moreover, rotating galax- 
ies tend to have smaller R^s than nonrotating galaxies 
(see Fig. [Tic). Therefore, our results also have conse- 
quences for the measured evolution of the tilt of the fun- 
damental plane and hence for the slope in the relation 
between M/L and either Af or a (as discussed previ- 
ously by van der Marel & van Dokkum 2007b). Evolution 
in the tilt of the f undamental plane has be e n observed 
by many autho rs ("van der Wei et al.' '2004'; 'Tre u et all 
2005 : Ivan der W cl et al. 2005; di Scrcgo Alighicr i et all 
2005t iJorgensen et al.l [2005). However, it is a mat- 



ter of debate to what extent sampl e selection effects 
contribute to the obse rved evolution (jvan der Wei et all 
l2005tl^eu et al.ir2005D . To illustrate the consequences of 
the r esults presented in t his pa per we revisit the analysis 
from Ivan der Wei et all ()200^ . They showed that the 
probability that the observed distribution of Afvir and 
Af/Lvir is drawn from a parent population with the same 
distribution as the local galaxy population, apart from 
a constant amount of luminosity evolution inferred from 
the most massive galaxies, is as low as 0.14%. We repeat 
the analysis but now using the distribution of AL joans and 
Af/Ljeans and find that this increases the probability to 
3.0%. In other words, the evidence for mass-dependent 
evolution of M/L remains strong, but it does become 
weaker. 

These arguments illustrate that the often observed evo- 
lution of the tilt of the fundamental plane is not necessar- 
ily entirely due to mass-dependent evolution of the M/L, 
i.e., downsizing. In order to determine to what extent the 
slope and the scatter of the relation between M/L and 
M evolve, still deeper observations are required to over- 
come the biases caused by sample selection effects in the 
surveys conducted over the past few years. In addition, 
it may be necessary to use dynamical modeling, as we do 
in this paper, to overcome the shortcomings of virial and 
fundamental plane mass estimates. 
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6. CONCLUSIONS 

We use the spatial information of our previously pub- 
lished VLT/F0RS2 absorption-line spectroscopy to mea- 
sure mean stellar velocity and velocity dispersion pro- 
files of 25 field early-type galaxies in the redshift range 
0.6 < z < 1.2, with median redshift z = 0.97. The 
kinematical profiles can be reliably measured even in the 
most distant galaxies. Rotation is detected in the ma- 
jority of the sample. Surface brightness profiles are de- 
termined from HST imaging. Two-integral solutions of 
the Jeans equations for oblate axisymmetric, constant 
axial-ratio models are calculated to interpret the data, 
taking into account line-of-sight projection, seeing con- 
volution, and pixel/slit binning. This yields for each 
galaxy the degree of rotational support, as quantified by 
the parameter fc, and the mass-to-lig ht ratio M/Ljcans 
(|van der Marel &: van Dokkumir2007af) . 

The rapidly rotating (fc > 0.6) and slow-rotating (k < 
0.6) galaxies in the distant sample tend to have different 
global properties. The rapidly rotating galaxies tend to 
have later morphological types (SO or Sa) and tend to 
be less luminous and more elongated. The slow-rotating 
galaxies tend to have earlier morphological types (E) and 
tend to be more luminous and rounder. The systematic 
variations with luminosity and axial ratio suggest that 
the correlation with morp hological type is rea l, in agree- 
ment with the findings of iMoran et al.l (|2007l ). Local E 
and SO galaxies show a more complex and overlapping 
set of kinematic properties (e.g., Cappellari et al. 2007). 
The contrasting clean separation of rotating and nonro- 
tating galaxies according morphological type observed at 
z ~ 1 is most likely partially artificial and must be re- 
lated to the fact that round/face-on SO galaxies are often 
misclassified as E galaxies, at high redshift even more so 
than at low redshift. 

We have compiled a local comparison sample to study 
evolution of the field early-type (E-I-SO) galaxy rotation 
rate. The distribution of k in our z ~ 1 sample is statis- 
tically indistinguishable from the distribution of {V/cr)* , 
a comparable measure of rotational support, in the local 
sample. The relative fraction of rotating galaxies does 
not change significantly between z ~ 1 (63% ± 11%) and 



the present (61% ± 5%). If rotation is taken to be a 
reliable indicator of morphological type, then this pro- 
vides evidence for an unchanging fraction of SO galaxies 
in the early- type field galaxy popu lation with r e dshift . 
This conflicts with the findings of IMoran et al.l (|2007l ) 
who did not correct for the effects of seeing convolution 
and pixel/slit binning on the measured rotation curves. 
It is possible that this may have led to an underestimate 
of the number of rapidly rotating field galaxies at z ~ 0.5 
in their study. 

We have compared the mass-to-light ratio M/Ljcans 
from our spatially resolved models to the values M/Lvir 
inferred by applying the virial theorem to globally av- 
eraged quantities. For E galaxies, which are generally 
slow-rotating, we find good agreement using a homol- 
ogy parameter /3 = 5, which is consistent with the 
value calibrated locally. So there is no evidence for non- 
homologous evolution of E galaxies out to z '^ 1. On the 
other hand, for elongated, rotating galaxies (which are 
often SOs) we find that M/Ljcans is on average ^ 40% 
higher than A//Lvir (see Figs. [Hand [7]). For local galax- 
ies this trend is not observed. This may hint at non- 
homologous evolution of this galaxy population (/3 in- 
creasing with redshift); it may suggest that our dynam- 
ical models produce biased results in these galaxies be- 
cause of the neglect of cold disks; or it may suggest that 
the virial formula for M/L produces biased results when 
applied to poorly resolved, rapidly rotating galaxies at 
high redshift. We cannot unambiguously identify the 
true cause, but to us the latter seems to be the most 
straightforward explanation. If so, then studies of fun- 
damental plane evolution overestimate both the amount 
of evolution and the evolution in the tilt of the fundamen- 
tal plane (and the relation between M and M/L), which 
is generally interpreted as evidence of down-sizing. 
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